Superconductivity and short range order in metallic glasses Fea;Nii_2;Zr2 
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In amorphous superconductors, superconducting and vortex pinning properties are strongly linked 
to the absence of long range order. Consequently, superconductivity and vortex phases can be stud- 
ied to probe the underlying microstructure and order of the material. This is done here from resis- 
tance and local magnetization measurements in the superconducting state of Fei,Nii_a;Zr2 metallic 
glasses with < a; < 0.6. Firstly, we present typical superconducting properties such as the critical 
temperature and fields and their dependence on Fe content in these alloys. Then, the observations 
of peculiar clockwise hysteresis loops, wide double-step transitions and large magnetization fluctua- 
tions in glasses containing a large amount of Fe are analyzed to reveal a change in short range order 
with Fe content. 



I. INTRODUCTION 



Since the synthesis of the first amorphous alloy super- 
conductors by vapor deposition by Biickel and Hilsch^ in 
1954, and the following, fabricated in the niid-1970s by 
electron beam evaporation by CoUver and Hammond", 
our understanding of superconducting phenomenon in 
this new class of superconductors has greatly evolved. 
Although these first attempts at making amorphous al- 
loys were plagued by the inconvenient instability of the 
amorphous phase at room temperature, they have served 
to expose the differences between superconductivity in 
amorphous materials and in their crystalline counterpart. 
Nowadays, several techniques based on the rapid cool- 
ing of the melt are used to fabricate stable amorphous 
alloys and the past 25 years have seen the publication 
of many studies about superconductivity in such metal- 
lic glasses, especially those composed of transition metal 
alloy s^'^'^'^'^'^. These studies have, among other things, 
discussed the importance of including effects due to spin 
fluctuations in the predictions of Tc, especially in alloys 
containing Ni, Co, or Fe. Additionally, they have iden- 
tified the important role played by microstructure and 
short range order (SRO) in amorphous materials. Since 
microstructure critically depends on the fabrication pro- 
cess and the melt cooling rate, values for different super- 
conducting characteristics such as the critical tempera- 
ture Tc from different labs for the same alloy composition 
often vary significantly. 

Vortices in type II superconductors can typically also 
reveal important information about material structure. 
In particular, defects and dislocations provide pinning 
sites which prevent vortex movement and enhance pin- 
ning. On the contrary, the absence of long range or- 
der in amorphous materials greatly decreases pinning 
properties. Moreover, several superconducting proper- 
ties, for instance the shape and width of the Bc2 transi- 
tion, or magnetization hysteresis loops are direct con- 
sequences of the homogeneity of the material and its 
fiux trapping capacities. Binary and pseudo-binary com- 
pounds composed of the early transition (ET) metal Zr 
and late transition (LT) metal Cu, Ni, Co and Fe in the 
form LT:rETi_:r and (LT°LT5_^)yETi_y have shown ex- 



cellent glass-forming abilities over a wide compositional 
range^iliSiiiSiii. In this article, we study such a pseudo- 
binary compound, namely Fea;Nii_a;Zr2 superconducting 
glasses with < x < 0.6, based on electric transport 
and local magnetization measurements. The high pu- 
rity and the amorphous nature of these alloys, conferring 
them extremely weak pinning properties, have previously 
allowed us to investigate transversely ordered dynamic 
vortex phasesi^ii^iiii. Here, we exploit the relationship 
between vortex pinning and material structure to reveal 
a change in SRO in this glass series. This conclusion 
is reached based on the observation of peculiar hysteresis 
loops and magnetization fluctuations in alloys containing 
a relatively large amount of Fe. Superconducting prop- 
erties of the metallic glasses and their dependence on the 
Fe content x are also presented. 



II. EXPERIMENTAL METHODS 

Alloy buttons are prepared by arc-melting appropri- 
ate amounts of the elemental constituents Fe (99.9%), 
Ni (99.999%), Zr (99.95%) under Ti-gettered argon at- 
mosphere in order to avoid oxidation. The buttons are 
re-meltcd 3 times to ensure homogeneity. Amorphous 
ribbons are then prepared by melt-spinning the alloy but- 
tons. Melt-spinning is performed in 40 kPa helium onto 
a copper wheel spinning at 50 m/s which ensures that 
the rapid cooling rate of 10^-10® K/s necessary for the 
formation of the amorphous phase is attained. The ab- 
sence of crystallinity was confirmed from the absence of 
constructive Bragg peaks in x-ray Cu Kq, diffraction (Fig. 
[1]). Indium contacts are soldered to the samples to per- 
mit electrical measurements in the standard four-probe 
technique. Resistance measurements are performed with 
a resistance bridge providing ac current at 15.9 Hz in a 
'''He refrigerator. The use of a dilution refrigerator was 
also required for measurements of the superconducting 
properties of Feo.6Nio.4Zr2 due to the low Tc below 0.3 
K. In both the '^He system and the dilution refrigerator 
a superconducting magnet with field capability up to 9 
T was used. The temperature is determined from cali- 
brated Cernox and RuO resistors. 
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FIG. 1: X-ray difFraction spectrum of Feo.3Nio.7Zr2 measured 
with Cu Kq radiation. The first peak is from diffraction from 
the glass support. The second peak is from diffraction from 
the sample. The position of this peak is determined from a 
Gaussian fit and the value is used in Ehrenfest's relation to 
obtain nnd. Inset: Near-neighbor distance for each alloy. The 
red line shows the mean nnd for the whole alloy series. 



A. Glass structure 

The question of the nature of the ordering in the amor- 
phous phase in the pseudo-binary Fe2:Nii_a;Zr2 series 
deserves particular attention. In the past, it has been 
tacitly assumed that the SRO characterizing the amor- 
phous structure of these alloys does not change upon 
substitution of Ni for Fe because these atoms have very 
similar sizes. Thereupon, various studies assuming con- 
stant glass structure in these, and similar glasses were 
undertakeni2iiiii^ii£ii2 to study the dependence of cer- 
tain effects on glass composition, independently of struc- 
tural change effects. Constant geometrical short range 
order (GSRO) across the series of alloys, mainly pro- 
vided by constant near-neighbor distances (nnd) , is read- 
ily verified from the position of the primary diffraction 
peak 9 in x-ray diffraction data and from the Ehrenfest 
relatioi*i^ii^ f — 0.6148A/ sin 6*, the mean near-neighbor 
distance f is evaluated, as shown in the inset of Fig.[l] us- 
ing A — 1.5405 A, the wavelength of Cu Kq radiation. As 
observed, the nnd in this alloy series is indeed constant, 
with a mean nnd of 2.998 A, which confirms constant 
GSRO in these alloys. This however does not necessarily 
imply that SRO is constant: since Fe and Ni have dif- 
ferent electronic structures, chemical short range order 
(CSRO), pertaining mainly to the atomic species of near- 
neighbors and their arrangement, cannot be assumed to 
remain constant. This question was previously investi- 
gated in these alloys using Mossbauer spectroscopy^^ but 
no change in SRO with x could be evidenced outside ex- 
perimental uncertainties. However, some results about 



superconductivity in this alloy series point to a transi- 
tion in SRO—; these will be shown and discussed later in 
this article. 



III. RESULTS 
A. Superconducting properties 

Initial interest in the study of superconductivity in 
the metallic glasses Fea;Nii_a;Zr2 was based on the as- 
sumed constant glass structure. Indeed, this would per- 
mit a study of the dependence of superconductivity on 
alloy composition, and more specifically the influence of 
spin fluctuations, independently of structure-dependent 
effectsiS,. Spin fluctuations, induced by the presence of 
magnetic atoms in the alloys, tend to demote supercon- 
ductivity by causing spin flips which break Cooper pairs. 
Consequently, one expects a suppression of the typi- 
cal parameters characterizing superconductivity, such as 
critical temperature T^ and upper critical field Bc2 with 
increasing Fe content. Such a behavior is witnessed in 
these alloys, as shown in Fig. [5^) and b). As can be seen, 
Tc of these alloys varies from 2.6 K to about 0.2 K with 
increasing Fe content. Tc is determined from resistance 
measurements in zero magnetic field and defined when 
the resistance reaches half the normal state value i.e. at 
R„/2. The values reported for the alloys are from several 
measurements on up to five different samples of each alloy 
composition. Among samples of the same alloy, typical 
Tc variations observed are smaller than 0.1 K; such a dis- 
tribution of TcS in an alloy is expected and inherent to 
the fabrication process. Indeed, the copper wheel used 
for melt-spinning becomes hotter in the process such that 
not the whole ribbon is cooled at exactly the same rate 
and the beginning of the ribbon can show significant dif- 
ferences in microstructure compared to the end of the 
ribbon; these differences are then reflected in supercon- 
ducting properties. Typically, annealing the samples at 
a temperature close to the glass transition temperature 
will remove these differences in microstructure along the 
length of the ribbon. However, doing this, we never- 
theless get a distribution of TcS in a single alloy. The 
augmentation of Tc from x=0 to 0.1 could result from an 
enhancement of the density of states at the Fermi level in 
Feo.1Nio.9Zr2 compared to NiZr2; this idea is supported 
by evidence from ultraviolet photoemission spectroscopy 
on binary alloys of Fe-Zr and Ni-Zr which have shown 
that the Fe d band lies closer to the Fermi level than the 
Ni d bandii . 

Similarly, the upper critical field Bc2 (Fig.[2]3)) was de- 
termined from resistance measurements performed as a 
function of magnetic field at a fixed temperature between 
0.33 and 0.35 K for the < a; < 0.5 alloys using a low 
driving current density J < 0.5 A/cm^. For the x = 0.6 
alloy, the temperature was below 0.1 K. The Bc2 transi- 
tion is defined at R„/2 and the values reported are from a 
few measurements on different samples. It will be shown 
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FIG. 2: a) Critical temperature and b) Upper critical field 
as a function of Fe content x in Fe2,Nii_2,Zr2. Error bars are 
statistical. For x=0.6, T < 100 mK during measurement of 

Bc2. 



later that the Bc2 transition exhibits large clockwise hys- 
teresis loops in the x = 0.5 and x = 0.6 alloys; in these 
cases, the value of Bc2 reported corresponds to the mean 
value obtained from up and down-going field sweeps, i.e. 
(5^2°""" - B^P)/2. It is found that 6^2 decreases from 
5.3 T to 1 T with increasing Fe content, and just like Tc, 
Bc2 increases for the alloy x=0.1 compared to x=0, again 
pointing to some promotion of superconductivity by the 
introduction of a small amount of Fc. 

In Table [H we report various physical and supercon- 
ducting parameters of the Fe2:Nii_3;Zr2 metallic glasses. 
The normal state resistivity p„ is calculated from resis- 
tance measurements performed at room temperature on 
long ribbons (> 30 cm) such as to minimize geometry- 
dependent effects, pn of the order 1.68 /zJlm is ob- 
tained for all alloys; these values are close to values 
reported for similar alloys^i^i^. Since pn is related to 
the GSRO, the constant pn throughout this composition 
range brings further confirmation of the constant GSRO. 
The slope of the upper critical field as a function of tem- 
perature close to Tc, '^^rf , is determined from a sin- 
gle set of measurements of the resistance as a function 



of temperature at different magnetic fields. We obtain 
^^§fr-\j, — —2.4 T/K as typical in amorphous alloys (see 

for instance^ '^'^1^^ ). We also report values for the pen- 
etration depth A(0), coherence length (0) (analogous 
to BCS ^0 but including correction due to short mean 
free path, as indicated by Gor'ko v^^'^^ ), and Ginzburg- 
Landau coherence length S^qi^ (0) and GL parameter k 
evaluated from expressions for superconductors in the 
dirty limi1>22. In NiZr2, we obtain A (0)= 0.87 //m and 
£^GL (0) ~ 8.1 nm; these values generally increase with 
Fe content to reach 3.04 fim and 16.2 nm respectively 
in x=0.6. This means that the vortex size and core in- 
crease substantially when going from x=0 to x=0.6. This 
should have important effects on vortex pinning proper- 
ties and vortex- vortex interactions. The GL parameter 
K is around 80 in these alloys, thus confirming that they 
are hard type-II superconductors. The mean free path is 
evaluated from 



/ = (3^2)1/3 



e^Pn 



2/3 



S_ 

Sf 



(1) 



where n,, is the free electron density and S/Sp is the ra- 
tio of the area of the free Fermi surface to that of a free 
electron gas of density rie. Both these quantities are es- 
timated as follows: is the ratio of the average number 
of electrons per atom outside closed shells to the atomic 
volume, i.e. ?ie = (f)^o^^- ^^^^ manner we have 
2.9 X 10^^ m"3 < rie < 3.2 x 10^^ m'^, with Ue decreas- 
ing with increasing Fe content. We also use S/Sp = 0.6 
as in Refi^; although if we were to use S/Sp = 1 as for 
a free electron Fermi surface the value of / would not 
change by an order of magnitude. As a result, we obtain 
I 2.8 Awhich is very close to the mean interatomic dis- 
tance (nnd = 2.998 A) in these amorphous alloys. As 
a result of the short mean free path, we obtain a dirti- 
ness parameter'^^ (0) /I above 20 which effectively con- 
firms that these amorphous alloys are in the dirty limit. 
We also evaluate the electron-phonon coupling parameter 
Xep — 0.6 using the McMillan equation 
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1.04 + ^* ln(ez3/1.45Tc) 
(1 - 0.62^*) In (615/1.457;) - 1.04 



(2) 



where we have used the Coulomb interaction parameter 
p* ~ 0.13 for polyvalent transition metala^^ and the De- 
bye temperature 8_d = 192.5 K as evaluated for NiZr2 ac- 
cording to Ref."^. This X^p makes the Fea;Nii_a;Zr2 metal- 
lic glasses in the weak to intermediate coupling regime. 



B. Structure inhomogeneity and CSRO 

Several physical properties of amorphous alloys de- 
pend on CSRO, for instance: the Curie temperature, the 
temperature coefficient of resistivity and superconduct- 
ing properties. In the case of superconductivity, this is 
because structural order directly influences the electron- 
phonon coupling parameter, resulting in a modification 
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TABLE I: Some measured and calculated physical and superconducting parameters. 



A 11 ^ 

Alloy 


r ^1 
p„ [fiil m] 


dB^o \rj-\ 1 ryl a 

df [^/-^^] 


A(0) [^mj 


(0) [nm] 


5gl (0) |nmj 


e 

n 




Bel [mL\-^ 


NiZr2 


1.68 ± 0.02 


-3.0 ± 0.2 


0.87 ± 5% 


6.7 ± 20% 


8.1 ± 4% 


96 ± 20% 


2.7 ± 40% 


0.175 ± 0.005 


Feo.1Nio.9Zr2 


1.68 ± 0.01 


-2.2 ± 0.4 


0.84 


7.5 


8.0 


101 


2.7 


0.197 


Feo.i5Nio.85Zr2 


1.62 ± 0.08 


-2.8 ± 0.1 


0.88 


7.2 


8.3 


106 


2.9 


N. A. 


Feo.2Nio.8Zr2 


1.69 ± 0.01 


-2.4 ± 0.4 


0.93 


8.0 


8.2 


99 


2.8 


0.295 


Feo.3Nio.7Zr2 


1.75 ± 0.02 


-2.8 ± 0.1 


0.98 


7.6 


8.8 


103 


2.7 


0.210 


Fco.33Nio.67Zr2 


1.84 ± 0.08 


-3.2 ± 0.1'' 


0.98 


7.0 


8.7 


86 


2.6 


N. A. 


Feo.36Nio.64Zr2 


1.72 ± 0.05 


-3.2 ± 0.1 


0.97 


7.1 


N. A. 


83 


2.7 


N. A. 


Feo.4Nio.6Zr2 


1.70 ± 0.02 


-2.6 ± 0.2 


1.02 


8.4 


9.6 


131 


2.8 


0.279 


Feo.5Nio.5Zr2 (1) 


1.69 ± 0.01 


-2.3 ± O.P 


1.22 


10.1 


12.3 


70 


2.8 


0.101 


Feo.6Nio.4Zr2 


1.67 ± 0.01 


N. A. 


3.04 


N. A. 


16.2 


N. A. 


2.9 


N. A. 



"The errors reported on | ^ consider the maximal and minimal slopes that could be obtained considering systematic errors on T 

and Bc2. 

''Obtained from A (0) = 1.05 X 10"^ iwY^^- 



"From (a (0) = 1.81 X 10" * 



dSc2 I 
dT 



-1/2 



'^From K = 3.54 X lO"* 



^From ^CL (0) = [^^fj^l^] ■ ^<^2 (0) is obtained from the extrapolation to T = of fits to the WHHM theorji2^ of Bc2 (T) data. 

-1I/2 

KBc2 
^ \ dT 

■'The % error on I is computed by considerering the effect on I of using a free electron-like Fermi surface ratio S/Sp = 1. 
f Measured at T ~ 0.35 K from local magnetization measurements^^. 

''Determined from resistance measurements as a function of magnetic field sweeps at different temperatures. 



of Tc for instance. Structural order can also influence su- 
perconductivity by acting on vortex-pin interactions, in 
particular defects, impurities and inhomogeneities, and 
can thus define the current-carrying capacities of the su- 
perconductor. For instance, it is the absence of long- 
range order in amorphous alloys which mainly determine 
their weak vortex pinning properties. Effects of struc- 
tural inhomogeneity are also commonly observed in su- 
perconducting properties, such as wide Bc2 or Tc transi- 
tions. In this section, we discuss some evidences of struc- 
tural inhomogeneity in the Fea;Nii_a;Zr2 with a large Fe 
content 0.4 < a; < 0.6 obtained from resistance and mag- 
netization measurements in the superconducting state. 
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FIG. 3: Width of the Bc2 transition for different Fea;Nii_a;Zr2 
alloys. In x=0.6, this includes both transitions 1 and 2. 



2. Clockwise hysteresis 



1. Width of transition 



A first obvious sign of the growth of structural inho- 
mogeneity with Fe content in these alloys is provided by 
an increase in the width of the Bc2 transition ABc2 = 
B(0.9R„)-B(0.1R„), as shown in Fig. [3l A double-step 
transition is even observed in some of the alloys x=0.5 
and 0.6 (Fig. 3]), an indisputable indication of inhomo- 
geneity. 



In type-II superconductors, the Bc2 transition between 
the normal state and the superconducting state often ex- 
hibits hysteresis due to Joule heating where more power is 
dissipated in the normal state than in the superconduct- 
ing state, and to flux pinning and trapping. The hystere- 
sis loop is then counterclockwise, i.e. the Bc2 transition is 
higher upon increasing the magnetic field than when de- 
creasing it. These types of hysteresis loops are observed 
here for alloys with < x < 0.4 (Fig. 2]) in which the Bc2 
transition is slightly lower upon decreasing (dotted lines) 
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FIG. 4: Resistance as a function of magnetic field for different 
alloys. The solid and dotted lines are for increasing and de- 
creasing B sweep respectively. The magnetic field was swept 
at a rate of 0.0147 T/s and the temperature was below 0.35 
K. 



the magnetic field than upon increasing (solid hnes) it. 
This can be contrasted to the large clockwise hysteresis 
loops seen at Bc2 in the alloys x = 0.5 and x = 0.6. All 
the X = 0.5 samples measured (7) show this wide clock- 
wise hysteresis loop, although only 2 of them show the 
double-step transition. The only x — 0.6 sample mea- 
sured shows a very broad Bc2 transition about 2 T wide 
including the two steps. A reversal of the hysteresis loop 
direction in x = 0.6 is observed between transition 1 and 
2 as identified in the figure; the uppermost transition 
exhibits the usual counterclockwise hysteresis loops. 

We have resistively measured the Bc2 transition in 
x=0.5 and 0.6 for diff'erent B field sweep rates. The re- 
sults, presented in Fig. [51 show an increase in the size of 
the hysteresis loops with increasing B sweep rate. This 
dependence on sweep rate provides evidence that a dy- 
namical process, such as vortex motion, is at the origin 
of the hysteresis loops. This is further supported by the 
data shown in the inset of Fig. [5l which shows the re- 
sistance as a function of time when the magnetic field 
sweep is paused in the middle of the Bc2 transition dur- 
ing an increasing and a decreasing magnetic field sweep. 
Time i = corresponds to the moment when the field 
sweep is paused. As can be seen, after t = the resis- 
tance initially keeps increasing (decreasing) over a short 
period of time for increasing (decreasing) B sweep, but 
it eventually reverses and decreases (increases) back to 
zero (a value close to what it was before the sweep was 
paused). Fitting these time dependences to an exponen- 
tial obtains a time constant r of 46.3 s and 43.8 s for the 
up and down-going field sweeps respectively, indicating a 
slow dynamical process. On the contrary to what is ob- 
served in magnetic field sweeps, no clockwise hysteresis 
loops are observed in temperature sweeps across Tc in a 
fixed magnetic field (data not shown). This brings fur- 
ther confirmation that the clockwise hysteresis loops are 
governed by a dynamical process involving vortex motion 
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FIG. 5: Width of clockwise hysteresis loop as a function 
of magnetic field sweep rate. AB is defined as the differ- 
ence in magnetic field when the resistance reaches 0.5Rn 
and up-going (B^j ) field sweep, i.e. 
) /2- Only the lowermost transition (tran- 
sition 2) is considered here in x=0.6. For Feo.5Nio.5Zr2, 
T=0.34 K and 1=30 For Feo.6Nio.4Zr2, T<0.1 K and 1=5 
/iA. The line is a guide for the eye. Inset: Resistance as a 
function of time when B field sweep is paused in the middle 
of the Bc2 transition. 1=0.03 mA, T=0.33 K. 



rather than some phase transition. 

In the literature, simulations of V-I characteristics 
in superconductors with inhomogeneous pinning poten- 
tials show clockwise hysteresis loop o^^i'^° . Such hysteresis 
loops result from dynamical effects and interplay between 
vortex trapping in the strong and weak pinning regions. 
As such, the size (or width) of the hysteresis loops is seen 
to depend on the sweep rate of the external variable with 
respect to which the loop is observed; in Refsi^ and^, 
this is the driving current and driving force. As is the 
case in our B field-induced hysteresis loops, the size of 
the loops increases in these Refsi^i^S for faster driving 
force sweep speeds. 

Based on these simulation results and our data, we pro- 
pose that the anomalous clockwise hysteresis loops ob- 
served at the Bc2 transition here arise due to the presence 
of inhomogeneities, some having stronger and weaker 
pinning properties, and thus resulting in an inhomoge- 
neous distribution of vortices in the superconductors. 
If we picture the sample as being composed of regions 
having stronger pinning properties than the surrounding 
medium, and if this medium also provides the connected 
path across the sample, we can explain the appearance 
of clockwise hysteresis loops in magnetic field sweeps and 
their absence in temperature sweeps as follows: As B is 
increased from 0, vortices will first penetrate in the main 
connected phase since it has the lowest pinning prop- 
erties and thus a lower energy barrier against flux en- 
try. Then, due to the elasticity of the vortex lattice, as 
the magnetic field is increased further it will be ener- 



6 



getically more favorable for the vortices to bend around 
the strong pinning regions'^^ and remain in the main con- 
nected phase. This will result in an inhomogeneous dis- 
tribution of magnetic flux in the superconductor with a 
larger flux density being present in the main connected 
phase, i.e. the phase of which the superconducting prop- 
erties are measured in resistance measurements. As a 
consequence, the Bc2 transition appears lower upon in- 
creasing the magnetic field than it would be if the fluxes 
were homogeneously distributed throughout the sample. 
A schematic representation of this process is presented 
in Fig. [B] On the contrary, upon decreasing the mag- 
netic field from above Bc2, fluxes tend to stay trapped 
in the strong-pinning regions but easily leave the weak- 
pinning main phase which again results in an inhomoge- 
neous distribution of vortices in the sample, but this time 
with the lowest vortex density in the main phase. This 
results in a resistively-measured Bc2 transition upon de- 
creasing the magnetic field that is higher than it should 
be. These two processes then result in clockwise hystere- 
sis loops. Similar phenomena have been observed previ- 
ously in granular supcrconductora^ii^i^^ and inhomoge- 
neous superconductors^i^SiM. 
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FIG. 6: Schematic representation of the resistive Bc2 tran- 
sition for up and down-going field sweeps with correspond- 
ing inhomogeneous vortex distribution in the superconductor. 
For increasing B field, more vortices pass through the weakly- 
pinned Fe-rich phase, avoiding the Ni-rich regions and the B"!" 
transition appears lower than it should be if the fluxes were 
homogeneously distributed in the whole sample (Wc2^). The 
inverse phenomena occurs upon decreasing B fleld, with vor- 
tices tending to remain in the more strongly pinned Ni-rich 
clusters such that Bff™" appears higher than W°i^ . 



This model explains why we witness an increase in 
the width of hysteresis loops with increasing B sweep 
rate: for slower sweep rate the vortices have more time 
to penetrate into the strong-pinning grains or to diffuse 
into the weak-pinning regions upon increasing and de- 
creasing the magnetic field before the measurement is 
taken. Therefore, after we apply a certain field, and 
by the time we take the measurement, the vortex dis- 
tribution has reached a more homogeneous configuration 



for slower sweep speeds and results in smaller hysteresis 
loops. A similar observation was made by Liu et al^ and 
Xu et at— from numerical simulations and experimen- 
tal measurements of fiux creep in a superconductor with 
inhomogeneous pinning properties. In these cases how- 
ever, the size of hysteresis loops was observed to increase 
with increasing driving current and driving force sweep 
rate, but the result is equivalent: changing the sweep 
rate amounts to changing the observation time window. 
Namely, for a slow sweep rate, our observation window 
is too late to observe the large inhomogeneity in the flux 
distribution. 

The results presented in the inset of Fig. [5] also lead 
to the conclusion that the vortex distribution is inhomo- 
geneous in these alloys. For instance, when the field is 
paused in an increasing magnetic field sweep, the resis- 
tance initially keeps increasing because fluxes easily and 
rapidly enter the weak-pinning phase; however it eventu- 
ally starts decreasing as no more fluxes are added (paused 
B field) and the fluxes in the weakly-pinned regions start 
to diffuse in the strong pinning grains, thus yielding a 
more homogeneous vortex distribution which brings the 
observed Bc2 transition closer to the real value. The op- 
posite takes place when the field is paused in a decreas- 
ing B sweep; the resistance change over approximately 
the same period of time is however smaller and indicates 
that the flux distribution is more homogeneous at Bc2 in 
a decreasing B sweep than in an increasing one. 

According to this model, no clockwise hysteresis loops 
are expected in temperature sweeps performed in a fixed 
external field, because in this case, the vortex density is 
fixed and its distribution across the sample remains the 
same as the temperature is swept up and down. 



3. Magnetization fluctuations 

As described in detail elsewher o^^i"^^ , local magneti- 
zation measurements were performed on these metallic 
glasses (0 < X < 0.5) using a 2-dimensional electron gas 
(2DEG) Hall probe. Fluctuations in the magnetization, 
increasing with Fe content (Fig. [7]), were observed and 
analyzed to reveal the presence of large vortex clusters 
of over 70 vortices in the superconductors with a large 
Fe content x > 0.4. It was also argued that the vortex 
bundles likely arise in these alloys because they are com- 
posed of two phases having different SRO, and thus dif- 
ferent pinning properties. In all alloys with a; > 0, both 
a Fe-rich and an Ni-rich phase exist, as evidenced from 
larger magnetization fiuctuations in x = 0.1 compared 
to X = 0, but the regions of the Fe-rich phase become 
larger and more numerous in x > 0.4; these results point 
to the existence of a structural phase transition close to 
X — 0.4, with alloys x < 0.4 having mostly NiZr2-like 
SRO and alloys with x > 0.4 having mostly FeZr2-like 
SRO. 

Summarizing the results presented in this section, 
many evidences indicate phase separation into Fe-rich 



7 



0.18 
0.15 
0.12 
0.09 
0.06 
0.03 
0.00 



-I- 



-I 



100 
80 

60 ^ 

40 

20 

-10 



0.0 



0.1 



0.2 0.3 
Fe content (x) 



0.4 



0.5 



FIG. 7: Black squares: Relative size of magnetization fluctu- 
ations for different glasses Fe2,Nii_2,Zr2. Red triangles: Num- 
ber of vortices in clusters related to the magnetization fluctu- 
ations. 



and Ni-rich regions in the Fej;Nii_2;Zr2 metallic glasses, 
particularly in the compositions with a large Fe content. 
Both the gradual increase in the width of the Bc2 tran- 
sition (Fig. [3]) and increasing size of magnetization fluc- 
tuations (Fig. [7|) with X point to an augmentation of the 
presence of inhomogeneities, most likely Fe-rich regions. 
The very broad ABc2, double-step transition, and ap- 
pearance of large clockwise hysteresis loops in a; = 0.5 
and 0.6 further supports this idea and provides indica- 
tion that Fe-rich regions have become critically large in 



these alloys, but some Ni-rich regions remain. 



IV. CONCLUSIONS 

In summary, we have presented the Fe content de- 
pendence of superconductivity in the metallic glasses 
Fea;Nii_a;Zr2. As expected due to the augmentation of 
spin fluctuations, Tc and Bc2 decrease with increasing x; 
this decrease becomes very pronounced in x=0.6. Im- 
portant progress was also made in our understanding of 
short range order in this amorphous alloy series. While 
we have shown that GSRO remains constant across the 
series as expected, numerous evidence, such as an in- 
crease in Bf,2 transition width with x, the appearance 
of double-step transitions and large clockwise hysteresis 
loops, and the observation of large fluctuations in mag- 
netization in high Fe containing alloys, point to the con- 
clusion that two phases having different CSRO exist in 
these alloys . According to these results, a Fe-rich phase 
having SRO resembling that of FeZr2 appears with a Fe 
content as low as x=0.1, but becomes critical around 
x=0.5 where important effects on superconductivity are 
witnessed. These results emphasize the extreme sensitiv- 
ity of superconductivity, and more specifically of vortices, 
as a probe of the microstructure of materials. 
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The dirtiness parameter is usually computed from the ra- 
tio of the BCS coherence length to the mean free path. 
However, we use the experimentally determined (0) be- 
cause it represents the real coherence length of our samples 
with consideration for the short mean free path in the dirty 
limit. 



